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ABSTRACT 
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Earlier classification analyses found three types of gamma-ray bursts (short, 
long and intermediate in duration) in the BATSE sample. Recent works have 
shown that these three groups are also present in the RHESSI and the BeppoSaX 
databases. The duration distribution analysis of the bursts observed by the Swift 
satellite also favors the three-component model. In this paper, we extend the 
analysis of the Swift data with spectral information. We show, using the spectral 
hardness and the duration simultaneously, that the maximum likelihood method 
favors the three-component against the two-component model. The likelihood 
also shows that a fourth component is not needed. 

Subject headings: gamma-rays: bursts, methods: statistical, data analysis 



Introduction 



Decades ago Mazets et al. (1981) and Norris et al. (1984) suggested that there might 
be a separation in the duration distribution of gamma-ray bursts (GRBs). Kouveliotou et 
al. (1993) found bimodality in the distribution of the logarithms of the durations. Today it 
is widely accepted that the phys ics of these two groups (short and long bursts — called also 
as Type I and Type II classes ( IZhang et al. 1 120071 : iKann et al. 1 120081 : IZhang et al. 1 12009 



Lii et al. 



2010lV) 



( iNorris et al. 



are different, and these two kinds o: 



2001; Balazs et al. 2003; Fox et al. 2005 



' GRBs are differe nt phenomena 
The angular 



sky distri bution of the short BAT SE's GRBs is anisotropic (jVavrek et al. Il2008[ ). In the Swift 



Kann et al. 



12008 



database ( ISakamoto et al. Il2008[ ). the measured re dshift distributions for the tw o groups are 



also different: for short bursts the me dian is 0.4 (jO'Shaughnessy et al. II2008I ) and for the 



long ones it is 2.4 (IBagoly et al. II2006I ) 



In the Third BATSE Cat alog flMeegan et a 



yses 
several papers 



Horvath 



20061: 



20021: 



Hakkila et al. 



Hakkila et al. 



12000 



19961) — using uni- and multi-variate anal- 



Horvathl ffl998[) and iMukheriee et al. I fll998h found a third type of GRBs. Later 



Balastegui et al. 



20031: Horvath et al. 



12001 



200 



Raianiemi fc Mahonen 



Borgonovo 112004 ; 



2002 



Horvath et al. 



Chattopadhyay et al. 1 120071 ) confirmed the existence of this third ("intermediate" 



in duration) group in the same database. The celestial distributi on of this third group 



Maghocchetti et al. II2003I : IVavrek et al. Il2008r ). 



in the BATSE sample is also anisotropic (|Meszaros et al. I l2000al Jbl: iLitvin et al. 1 12001 



Recent works analy zed the Swift (IHorvath et a 



(IRlpa et al. I l2008l 120091 ) and BeppoSaX flHorvath 



20081 iHuia et al. 



mm, RHESSI 
20091 ) data, respectively. They have 



found the intermediate class in all the three satellites' data: in the Swift database the one- 



3 



dimensional maximum likelihood (ML) analy sis of the durations has proven the existence of 



the BeppoSAX database ( iFrontera et al. 



these three subgroups (IHorvath et al. Il2008[): a preliminary study by the same method o f 



20091 ) gave support for this class (IHorvath II2009I ): 



in the RHESSI database two methods led to the same results - the same one-dimensional 
ML method of the durat ions and the bivariate ML method using both duration and hardness 
feipa et al. Iboosi boool ). 



A method to infer t he physical origin o f GRBs was developed recently (jZhang et al. 



2007: 


Kann et al. 


2008: 


Zhang et al. 


2009) 



ration are used as the differentiation criteria. Such a scheme only result in two major types 
of GRBs (Type I and Type II). On t he other hand, G RBs may be classified using different 



parameters other than duration (e.g., 
discussion section. 



Lii et al. 



We shall discuss these more in the 



taneously; 



Ripa et al. 



Horvath et al. I (120061) an alyzed the BATSE data using duration and hardness simul- 



(120091 ) studied the RHESSI data with t he same configu r ation . The 
bivariate analysis on the Swift database has not been done yet. IHorvath et al. I (120081 ) only 
provided the one-dimensional analysis of the durations. 

Hence, to get a complete picture, one has to analyze the Swift data also - using both 
the duration and hardness simultaneously - with the bivariate ML method. This is the aim 
of this paper. 

The paper is organized as follows. Section 2 briefly summarizes the method of the 
two-dimensional fits. Section 3 defines the sample. Section 4 deals with these fits in the two- 
dimensional parameter space and confirms the reality of the intermediate group. Section 5 
discusses the physical differences between the classes. Section 6 contains the discussion and 
Section 7 summarizes the conclusions of this paper. 



2. The mathematics of the method 



When studying a GRB distribution, one can assume that the observed probability dis- 
tribution in the parameter space is a superposition of the distributions characterizing the 
different types of bursts present in the sample. Using the notations x and y for the variables 
(in a two-dimensional space) and using the law of full probabilities (IRenyi Ill962l ). one can 
write 



1=1 
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In this equation, p{x, y\l) is the conditional probabihty density assuming that a burst belongs 

k 

to the Ith class, pi is the probability for this class in the observed sample = 1); where 

1=1 

k is the number of classes. In order to decompose the observed probability distribution 
p{x,y) into the superposition of different classes, we need the functional form of p{x,y\l). 
The probability distribution of the logarithm of durations can be w ell fitted by Ga ussian 
distributions, if we restrict ourselves to the short and long GRBs (IHorvath I Il998[ ). We 
assume the same also for the y coordinate. With this assumption we obtain, for a certain 
l-th class of GRBs, 



exp 



p{x,y\l) 



2(1 -r^ 



1 



2-Kax(yyVl — 
{x-a^y , {y- 



X 



2r{x - ax){y - ay) 



a^Oy 



(2) 



where the meaii s, a^, Cy are the dispersions, and r is the correlation coefficient 

(iTrumpler Sz Weaver Hence, a certain class is defined by five independent parameters, 

and r, which are different for different /. If we have k classes, then we have 



Cly, (Tj 



0", 



(6A;— 1) independent parameters (constants), because any class is given by the five parameters 

k 

of Eq.(l2]) and the weight pi of the class. One weight is not independent, because YIpi ~ ^■ 

1=1 

The sum of k functions defined by Eq.([2]) gives the theoretical function of the fit. 

By decomposing p{x,y) into the superposition of p{x,y\l) conditional probabilities one 
divides the original population of GRBs into k groups. Decomposing the left-hand side 
of Eq.([T]) into the sum of the right-hand side, one needs the functional form of p{x,y\l) 
distributions, and also k has to be fixed. Because we assume that the functional form is a 
bivariate Gaussian distribution (see Eq.([2])), our task is reduced to evaluating its parameters, 
k and pi. 



Balazs et al. I (|2003[ ) used this method for k = 2, and gave a more detailed description of 



the proced ure. However, that pap er used fiuence instead of hardness, and used the BATSE 
data. Also iHorvath et al. I (|2006[ ) used this method analyzing the BATSE data. Here we 
will make similar calculations for k = 2, k = 3 and k = A using Swift observations. 



3. The sample 



In the Swift BAT Catalog fISakamoto et al. 1120081 ) there are 237 GRBs, of which 222 have 
duration information. Following the same procedure of data reduction we have extended this 
sample with all the bursts detected until mid 2008 December (ending with GRB 081211a). 
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Our total sample thus comprises the first four years of the Swift satellite (since the dete ction 



of its first burst GRB 041217) and includes 325 bursts. 222 from lSakamoto et al. I fl2008[ ) and 
103 reduced by us. The data reduction was done by using HEAsoft v. 6. 3. 2 and calibration 
database v. 20070924. For light curves and spectra we ran the batgrbproduct pipeline 0We 
fitted the spectra integrated for the duration of t he burst with a power-l aw model and a 
power-law model with an exponential cutoff. As in ISakamoto et al. I (120081 ) we have chosen 
the cutoff power-law model if the of the fit improved by more than 6. 

For calculating the hardness ratio, we have chosen fiuence 2 (25 — 50keV) and fiuence 
3 (50 — lOOkeV) and the hardness is defined by the HR = F3/F2 ratio. Figure [1] shows the 
distribution of GRBs in the {IgTgo; Ig-ff-R} plane, where the fits were made for x = IgTgo 
and y = Ig HR. 



4. Bivariate ML fitting 

In order to find the unknown co nstants in Eq.(f2l), we use the maximum likelihood (ML) 



procedure of parameter estimation (IBalazs et al. 



20031 ). Assuming a set of N observed 



= 1, . . . , A^) values (A^ is the number of GRBs in the sample for our case, which 
here is 325) we can define the likelihood function in the usual way, after fixing the value of 
k, in the form 



TV 



^lnp( 



(3) 



where p(xi, yi) has the form given by Eq.([T]). Similarly to what was done by iBalazs et al. 



( 120031 ) and iHorvath et al. I ( 120061 ). the EM (Expectation and Maximization) algorithm is 



used to obtain the ax, ay, ax, (Jy, r, and pi parameters at which L reaches its maximum value. 
We made the calculations for different values of k in order to see the improvement of L, as 
we increase the number of parameters to be estimated. Tables [T]|2] summarize the results of 
the fits for k = 2 and k = 3. Figures 2 and 3 show the results in the {IgTgo; ^SHR} plane. 



^ |http: / /heasarc. nasa.gov/lheasoft/ 



Table 1: Results of the EM algorithm in the {IgTgo; ^gHR} plane, k = 2 



I 


Pi 


ax 


ay 


0"x 


(Ty 


r 


1 


0.082 


-0.383 


0.256 


0.602 


0.114 


0.071 


2 


0.918 


1.628 


0.096 


0.516 


0.117 


0.226 
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Fig. 1.— Distribution of = 325 GRBs in the {IgTgo; IgHR} plane. 
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Fig. 2. — Distribution of = 325 GRBs in the {IgTgo; ^sHR} plane. The different symbols 
mark bursts belonging to the short (red +) and long (green x) classes, respectively. 
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The confidence interval of tlie estimated parameters can be given on tlie basis of tlie 
following theorem. Denoting by L^ax and Lq the values of the likelihood function at the 
maximum and at the true va lue of the parameters, re spectively, one can write asymptotically 



as the sample size N ^ oo ( iKendall &: Stuart Ill983[ ) 



2{L„,ax - Lo) ^ xl (4) 
where m is the number of estima ted parameters (m = 6k — 1 in our case), and Xm ^^e 



usual m-dimensional function (ITrumpler fc Weaver Ill953l ). Moving from k = 2 to k = 3 



the number of parameters m increases by 6 (from 11 to 17) and Lmax grows from 506.6 to 
531.4. Since Xi7 = Xu + Xe increase in L^ax by a value of 25 corresponds to a value of 
50 for a Xq distribution. The probability for xl > 50 is very low (10^^), so we may conclude 
that the inclusion of a third class into the fitting procedure is well justified by a very high 
level of significance. 

Moving from = 3 to = 4, however, the improvement in L^ax is 3.4 (from 531.4 to 
534.8) corresponding to Xq ^ 6.8, which can happen by chance with a probability of 33.9 %. 
Hence, the inclusion of the fourth class is not justified. We may conclude from this analysis 
that the superposition of three Gaussian bivariate distributions - and only these three ones 
- can describe the observed distribution. 

This means that the 17 independent constants for A; = 3 in Table [2] define the parameters 
of the three groups. We see that the mean hardness of the interr nediate class is v e ry low - the 



third class is the softest one. This is in a good agreement with iHorvath et al. I (|2006[ ). who 
found that the intermediate duration class is the softest in the BATSE database. In that 
database, 11% of all GRBs belonged to this group. In our analysis, p2 = 0.296; therefore, 
30% of the Swift bursts belong to the third group. 



5. Separation of GRBs into the classes 

Based on the calculations in the previous paragraph, we resolved the p{x, y) probability 
density of the observed quantities into a superposition of three Gaussian distributions. Using 



Table 2: Results of the EM algorithm, k = 3 



I 


Pi 


dx 


ay 


(Tx 


ay 


r 


1 


0.079 


-0.426 


0.259 


0.576 


0.114 


0.120 


2 


0.296 


1.076 


0.025 


0.376 


0.129 


-0.004 


3 


0.626 


1.882 


0.130 


0.350 


0.093 


-0.237 
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Fig. 3. — Distribution of = 325 GRBs in the {IgTgo; Ig HR} plane. The different symbols 
mark bursts belonging to the short (red +), intermediate (blue *) and long (green x) classes, 
respectively. 
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this decomposition, we can cl assify any observed GR B into the classes represented by these 



groups (this is similar to the iHorvath et al. I (120061 ) work dealing with the BATSE data 



In other words, we develop a method allowing us to obtain, for any given GRB, its three 
membership probabilities, which define the likelihood of the GRB to belong to the short, 
intermediate, and long groups. The sum of these three probabilities is unity. For this purpose 
we define the following Ii{x,y) indicator function, which assigns to each observed burst a 
membership probability in a given / class as follows: 

p{x,y\l)pi 

h{x,y) = - . (5) 

1=1 

According to Eq.([5]), each burst may belong to any of the classes with a certain proba- 
bility. In this sense, one cannot assign a given burst to a given class with absolute certainty, 
but with a given probability. This type of classification is called a "fuzzy" classification 



( iMcLachlan fc Basford 1119881 ). Although, any burst with a given [x,y] could be assigned 
to all classes with a certain probability, one can select that / at which the Ii{x,y) indicator 
function reaches its maximum value. For k = 3, Figure 3 shows the bursts' distribution in 
the {IgTgo; IgiY-R} plane indicating the group memberships. For the community, we make 
the list of the membership probabilities available on the internet. |§ Table [3] contains the six 
GRBs identified as intermediate having redshift information. 

One can demonstrate the robustness of classification by comparing the results obtained 
from the {IgTgo; Ig-ff-R} and the {lgT^o;\gHR} planes. A cross tabulation between these 
two classifications is given in Table |H According to the table, the short and long classes 
correspond within a few percent to the respective groups obtained from the other classifica- 
tion. We even have three zeros, which means, for example, no short burst, classified with 



The list of the membership probabihties can be found at 



http://itl7.elte.hu/^veresp/sak325T90H32CL3.res, 



Table 3: The six intermediate GRBs identified with redshift. 



GRB id. 


3rd type prob. 


redshift 


050525 


0.78 


0.606 


050922C 


0.71 


2.17 


060206 


0.88 


4.048 


071117 


0.82 


1.331 


080913 


0.67 


6.695 


081007 


1.00 


0.5295 
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{IgTgo; Ig-ff-R}, was identified as long with the {IgTso; Igi^-R} classification and since the 
other number is also zero, this applies vice versa. Consequently, the robustness of the short 
and long groups is well established. On the other hand, the population of the intermediate 
group is less numerous when classifying in the {IgTso; Ig HR} plane than in the other one. 
Table m clearly shows that the {lgTgo;lg HR} classification, except for three (two for long 
and one for short), contains all GRBs assigned to the intermediate group by {IgTso; Ig -ff-R}. 

The intermediate bursts classified in the {IgTgo; Igi^-R} plane indicated GRBs from 
the short and 24 from the long group identified in the other plane. This moderately high 
number of indicated bursts clearly shows that a slight variation of the parameters of the 
Gaussian distribution representing the intermediate group results in a moderate change in 
the number of identified objects in this group. Comparing the number of GRBs belonging 
to the intermediate group, one gets 86 and 65; 86 are identified as intermediate with the 
Tgo analysis and 65 are identified as intermediate with the T^o analysis (see the numbers in 
Table SD. 

If one assigned the burst to that group that had the maximum membership probability, 
a slight change in the parameters of the corresponding Gaussian distribution may move 
the GRB to an other group. On the contrary, the fuzzy classification assigns membership 
probability to all of the bursts. Hence, a small variation of the parameter gives a small 
variation in the estimated number of bursts in the intermediate group obtained by summing 
the membership probabilities of all GRBs in the sample. There is a further issue to be 
considered here. The duration (Tgo) is dependent on the energy range at which we are 
measuring. This results in a further fuzziness factor in the determination of the membership. 



Table 4: The cross-tabulation of the two classifications (with Tgo — HR and T50 — HR). 



T50 \ T90 


short (Tgo - HR) 


interm. (Tgo — HR) 


long (Tgo - HR) 


total 


short (T50 - HR) 


24 








24 


interm. (T50 — HR) 


1 


62 


2 


65 


long (T50 - HR) 





24 


212 


236 


total 


25 


86 


214 


325 
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Discussion 



6.1. Redshift Distributions 



The cumulative redsh ift distributio n of t he three populations is shown in Figure HJ 
Redshifts were taken from lAmati et al. I (|2008[ ). Only a subset of the classified bursts had 
redshift information and we considered bursts where the probability of belonging to a given 
population is higher than 97%. This means 6 short, 9 intermediate, and 50 long GRBs. 
The long and short population redshift distributions are significantly different (99.4% sig- 
nificance). The intermediate GRBs redshift distribution is clearly between the short and 
long redshift distributions, which could mean that they are further than the short bursts 
and closer than the long ones. However, probably owing to the small number of data points 
the difference is not significant. We have tried several statistical tests, such as Kolmogorov- 
Smirnov, Wilcoxon, and Mann- Whitney. None of them showed high significance; the best 
one was 92%. Therefore, we are not able to prove that the long and the intermediate bursts 
came from different redshift distributions. A more extensive analysis on these aspects as 
well as the af terglow properties of the interm ediate group of events will be presented in a 
parallel work (jde Ugarte Postigo et al. Il2010h . 



6.2. The -Epcak — ^iso (Amati) Relation 



Once we classify the bursts, it is also possible to investiga te their properties in the 
context of the -Epeak — E^so or Amati-relat ion (lAmati et al. 1120081) in the case of bursts with 
measured redshift. The Amati-relation (lAmati et al. 1 120021 ) is a correlation between the 
rest-frame peak energy of the GRB spectrum (-Epcak) and the isotropic-equivalent energy 
release of the burst (E'iso). 

Again, only a fraction of the populations had redshift and could be placed on the 
Epeak — Eiso plane. We have 18 from the long population and 6 from the intermediate. Both 
groups seem to follow the same relationship. As the Amati-relation is not valid for the short 
population, the intermediate bursts are more closely related to the long population than to 
the short class. 

Intermediate bursts do not populate the most energetic regime of the -Epeak ~ -^'iso plane 
unlike the long bursts (see Figure [3]). They tend to have lower isotropic energies compared 
to the long population. The small number of data points makes hard to give firm assertions 
at this time. Also, there is no significant clustering of intermediate bursts on this plane. 
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Fig. 4. — Cumulative redshift distribution of the three classes. The continuous line is the 
short, the dashed line is the intermediate, and the dotted line is the long population. 
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Fig. 5. — Epcak — Eiso relation of the classified bursts. Red symbols represent long GRBs, 
and intermediate bursts are plotted in blue. The y-axis is the 10-base logarithm of ii^peak in 
keV, and the x-axis is the 10-base logarithm of Eiso in units of 10^^ erg/s. 
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6.3. Relation to Other Classifications 



One ca n compare this classi fication to the Type I/II method put forward by lZhang et al 



(2007 


) and 


Zhang et al. 


(2009) 



The Type 1/11 scheme looks for signatures of the binary 
merger and the collapsar scenario (as association with supernova, host galaxy properties, 
spectral lag etc.) in bursts and classifies them accordingly. Currently these two scenarios 
are thought to be the most probable progenitors of GRBs. This method uses a wide range 
of observations for classification (some deemed more d ecisive than others) and as such only 



a fraction of the bursts can be assigned Type 1 or 11. IZhang et al. I (120091 ) publish a table 
with the most certain members called the Gold sample. Though this scheme allows only for 
two classes it is worth checking their membership in our classification scheme. 

There are five bursts of Type 1 in the Gold sample. GRB050709 was detected by HETE 
11 so it is not included in our sample. The other four GRBs are listed in Table [5l This table 
c ontains the group membership probabilities of the bursts from the Type 1 Gold sample of 



Zhang et al. I (|2009[ ). which clearly shows none of them belong to the intermediate class. 



The lightcurves of the last three of the bursts in Table [5] have a short-hard spike followed 
by an extended emission tail. If we disregard the extended emission we get a different 
duration. The intermediate group membership probability increases in this case to ~ 0.95, 
0.98 and ~ 0.45 for 050724, 060614 and 061006 respectively. 061006 has a ~ 0.55 probability 
of belonging to the short group. Care must be taken with these probabilities as only a whole 
new classification with the new durations would yield the correct membership probabilities. 

GRB 060218 and X-ray fiash 080109 have associated supernovae. Unfortunately, there 
is no durati on available for an y of them. For 060218, a minimum duration of ~ 2700 s can be 
established (jSakamoto 1 120061 ). and the logarithm of hardness ratio of IgHR = —0.15 can be 
derived. These values place this burst to the extreme right of the hardness-duration diagram 
and thus it belongs to the long population (though it is very soft). 080109 was detected with 
XRT and BAT only measured upper limits; therefore, we cannot infer anything about its 



Table 5: Membership probabilities of the Type 1 Gold sample of IZhang et al. I (|2009[ ) 
(columns 2-4). The memberships calculated without the extended emission (EE) are in 
the last three columns (marked with a *^ 



GRB name 


P( short) 


P(interm.) 


P(long) 


P* (short) 


P*(interm) 


P*(long) 


050509B 


1 








1 








050724 





0.02 


0.98 


0.02 


0.95 


0.03 


060614 





0.02 


0.98 


0.01 


0.98 


0.01 


061006 





0.01 


0.99 


0.53 


0.44 


0.03 
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group membership. 



Conclusion 



In the BATSE analysis iHorvath et al. I (l2006l ) found the intermediate type of burst to be 
the softest. To be sure that the classifications are free of instrumental effects, it is important 
to compare the results obtained with datasets from different satellites. This is one of the 
main aims of this paper. We find that the bursts observed by Swift can be divided (in 
the duration hardness plane) by three groups and only three groups as happened with the 
BATSE sample. 



Our results, summarized in Table 2, are very similar to iHorvath et al. I (120061 ) results 
obtained using the BATSE data. This indicates that the two satellites are observing a similar 
population of bursts. The relative size of the different groups will differ with the detector 
parameters: the Swift observations are more sensitive for soft and weak bursts; hence, the 
observation probability of the intermediate group members is slightly enhanced compared to 
BATSE (from 11% to 30% ). 

An important question that must be answered in this context is whether the intermediate 
group of GRBs, obtained in the previous paragraph from the mathematical phenomenological 
classification, really represents a third type of burst physically different from both the short 
and the long population. 

To infer the physical origin of a GRB, one must collect more direct information about 
the GRB progenitors. For e xample, most long GR Bs are found to have irregular host galaxies 
with intense s tar fo rmation (jPruchter et al. Il2006[ ) and some are associated with a supernova 
(jPian et al. I ( 120061 ) . and references th erein). Some short GRBs are associated with nearby 
galaxies with low star fo rmation rate ( iGehrels et al. I l2005t iFox et al. I l2005l : iBerger et al. 
20051 ; IZhang et al. 1 120071 ). which point toward a possible origin of compact star mergers. 



Since one has only two main types of suggested progenitors, massive star progenitors 
for long (Type II) and compact star mergers for short (Type I), there are two possibilities: 

A, There are only two physically different types of GRBs, and the statistically significant 
third group belongs to one of them. 

B, The third group is physically real, and one should look more carefully at the new 
observations to find this new type of progenitor. 



According to this paper's analysis, the intermediate GRBs are the softest among the 



- 16 - 



three classes. This different small mean hardness and also the different average duration 
suggest that the intermediate group should also be a different phenomenon, that is, both in 
hardness and in duration the third group differs from the other two. 

To sum up, our bivariate ML method confirmed the existence of the third intermediate 
subgroup on the high significance level in the Swift database. Existence of other groups 
is not supported. It is conjectured that the GRBs of intermediate subgroup are physically 
different phenomena. 

This research was supported in part through OTKA K077795 grant, by the GAUK 
grant No. 46307, by the GACR grant No. P209/10/0734, by OTKA/NKTH A08-77719 and 
A08-77815 (ZB), by the Research Program MSM0021620860 of the Ministry of Education of 
the Czech RepubUc (AM), by ASI grant SWIFT 1/011/07/0, by the Ministry of University 
and Research of Italy (PRIN MIUR 2007TNYZXL) (AdUP), by a Bolyai Scholarship (IH), 
and by an ESO fellowship (AdUP) . Thanks are due to the valuable remarks to Jakub Ri'pa. 
The detailed and constructive suggestions and remarks of the anonymous referee significantly 
helped to improve the scientific merit of this paper. 
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